x <- c(0.6)
y <- c(0.2)
alpha <- 1E-3
beta <- 1.2E-2
gamma <-7E-3
delta <- 2E-2


for (i in 1:4000) {
  lastx <- x[length(x)]
  lasty <- y[length(y)]
  dx = lastx*(alpha - beta*lasty)
  dy = -lasty*(gamma - delta*lastx)
  x <- c(x, lastx + dx)
  y <- c(y, lasty + dy)
}

plot(x, type="l", col="green", ylim=c(0,max(x)))
lines(y, type="l", col="black")

#plot(x,y, type="l", col="gray", xlab="dienas",ylab="zvēri")
